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Abstract 

The non-mesonic weak decay of polarized A-hypernuclei is studied with a microscopic diagram- 
matic formahsm in which one- and two-nucleon induced decay mechanisms, AA^ — )• A^A^ and 
ANN — > NNN, are considered together with (and on the same ground of) nucleon final state 



I interactions. We adopt a nuclear matter formalism extended to finite nuclei via the local density 

' approximation. Our approach adopts different one-meson-exchange weak transition potentials, 

. while the strong interaction effects are accounted for by a Bonn nucleon-nucleon interaction. We 

also consider the two-pion-exchange effect in the weak transition potential. Both the two-nucleon 
induced decay mechanism and the final state interactions reduce the magnitude of the asymmetry. 
The quantum interference terms considered in the present microscopic approach give rise to an 
opposite behavior of the asymmetry with increasing energy cuts to that observed in models de- 
scribing the nucleon final state interactions semi-classically via the intranuclear cascade code. Our 
results for the asymmetry parameter in j^C obtained with different potential models are consistent 
with the asymmetry measured at KEK. 

PACS numbers: 21.80.+a, 25.80.Pw. 
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I. INTRODUCTION 



The study of hypernuclear physics provides the main source of information on the baryon- 
baryon strangeness-changing weak interactions. In particular, the non-mesonic weak decay 
of A-hypernuclei has shown two challenging issues presenting some puzzhng character [1, 2]. 
First, we must mention the disagreement between theory and experiment for the ratio Vn/Vp 
between the rates for the An nn and the Kp — )■ np non-mesonic weak decay processes. 
Another more recent problem concerns the asymmetry in the proton emission from the 
non-mesonic weak decay of polarized hypernuclei, which is our main concern in the present 
contribution. 

A A-hypernucleus can be produced with some degree of polarization. Indeed, the 
n(7r"'", K^)k reaction has been used [3] at p^r = 1-05 GeV and small K'^ laboratory scatter- 
ing angles to produce hypernuclear states with a substantial amount of spin-polarization, 
preferentially aligned along the line normal to the reaction plane which identifies the po- 
larization axis. The dominant decay mechanisms for a polarized A-hypernucleus are the 
following neutron- and proton-induced processes: 



It turns out that the number of protons emitted parallel to the polarization axis is different 
from the same quantity measured in the opposite direction. This asymmetric proton emission 
is a consequence of the interference between the parity-conserving and the parity-violating 



Let us denote with a\ the intrinsic asymmetry, arising from the one-nucleon induced 
(lA^) decay in Eq. (2). The one-nucleon induced decays take place within the nuclear 
environment and the resulting nucleon pairs can interact strongly with others nucleons of 
the medium before any nucleon leaves the nucleus and is detected. As a result of these 
final state interactions (FSI), the asymmetry measured in an experiment, aj^f, differs from 
the intrinsic value, a]^ . Most of the theoretical models result in a negative and rather 
mass-independent intrinsic asymmetry. Instead, data favor a small ajjf, compatible with a 
vanishing value, for both ^He and j^^C. This shows a clear disagreement between aj^^ and 



An — > nn 



Kp — >■ np. 



(1) 
(2) 



terms in the Kp — > np weak transition potential [4]. 
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The reason for this disagreement can be twofold. It can be originated by the weak 
decay mechanism itself, which might require some improvement and the consideration of 
additional two-nucleon induced processes, and it may also be due to nucleon FSI. Let us 
start by considering the various mechanisms which contribute in the evaluation of a]^ . The 
theoretical models based on one-meson-exchange potentials (OME) [5-10] and/or direct 
quark mechanisms [11] predict values in the range from —0.7 up to —0.4. By using an 
effective field theory approach, a dominating central, spin- and isospin-independent contact 
term was predicted in [12] which allowed the authors to reproduce the experimental total 
and partial non-mesonic decay widths for ^Hc, and ]^C, and the asymmetry parameter 
for ^He. Motivated by this work, a scalar-isoscalar cr-meson-exchange was added to a 
(tt + i^)-exchange weak model also including a direct quark mechanism [13]. Similarly, the 
(T-meson was considered together with a full OME weak potential in [14]. Although the 
addition of the cr-meson may improve the calculation of a]^ , it turned out to be not enough 
to reproduce consistently all the decay data despite the freedom introduced by the unknown 
coupling constants of the cr-meson. Later, the OME weak potential was supplemented by 
the exchange of (uncorrelated and correlated) two-pion pairs [15]. The two-pion-exchange 
potential was obtained from a chiral unitary approach in a study of the nucleon-nucleon 
interaction [16] and was adapted to the weak sector in [17] for a study of the non-mesonic 
decay rates, while the calculation of the asymmetry was also carried out in [15]. The two- 
pion exchange mechanism turned out to introduce a significant central, spin- and isospin- 
independent AN — > nN amplitudes and gave rise to a good reproduction of the entire set 
of decay rates and asymmetry data for ^He and ]^C. 

We now briefly comment on the effect of FSI on the asymmetry parameter. First of 
all, it should be noted that from a strictly quantum-mechanical point of view the only 
observables in non-mesonic weak decay are the total non-mesonic decay width, Fnm, the 
spectra of the emitted nucleons and the asymmetry [18]. It is the action of FSI which 
prevents the measurement of any of the non-mesonic partial decay rates and of the intrinsic 
asymmetry aj^ . The link between theory and experiment for both F^/Fp and a]^ , is not 
straightforward, since it is strongly dependent on FSI. For instance, to obtain the F„/Fp 
ratio from experiments, one should proceed to a deconvolution of the nucleon rescattering 
effects contained in the measured nucleon spectra [19], which requires the use of a theoretical 
approach for FSI. For the asymmetry parameter the situation may seem more direct, as 
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experimental data for ajf are available. However, for a direct comparison with experiment, 
one must calculate the asymmetry ajf , which also requires the inclusion of FSI effects. Only 
a couple of approaches [9, 15] calculated this observable in an appropriate way. However, 
both these calculations adopted an hybrid approach consisting in a shell model for describing 
the weak decay and a semi-classical intranuclear cascade (INC) model, for simulating FSI. 
The only kind of FSI effects considered within the finite nucleus approach of [9, 15] are 
those between the two nucleons emitted in the non-mesonic decay, which are represented 
by a wave function describing their relative motion under the influence of a suitable NN 
interaction [7] . Although a discrepancy with data still remains for proton emission spectra 
[20, 21], one can safely assert that a formalism which takes care of FSI leads to a good 
agreement between theory and experiment concerning F^/Fp and a^. 

In the present contribution we evaluate the asymmetry ajf employing an alternative 
approach to the hybrid one of [9, 15]. In our microscopic diagrammatic approach, which 
was developed in [20, 22, 23], both the weak decay and the nucleon FSI are part of the same 
quantum-mechanical problem and are thus described in a unified way. Therefore, the present 
formalism has a self-consistency that is not present in previous approaches. The calculation 
is first performed in nuclear matter and then extended to finite hypcrnuclei by means of 
a local density approximation. Clearly, our nuclear matter wave functions are less realistic 
than the shell model ones. However, FSI are relevant and our quantum-mechanical approach 
describes them more rehably than the INC. In [20] we showed that quantum interference 
terms in the FSI are very important in the calculations of the observable spectra for the 
emitted nucleons. Moreover, in the same work we have called attention on the fact that 
pure (i.e., non-quantum interference terms) FSI terms and two-nucleon induced (2A^) decay 
contributions originates from two different time-orderings of the same Feynman diagrams 
at second order in the weak transition potential. 

Another contribution of the present work is the consideration for the first time of the 2A^ 

— * 

decays, ANN — >■ nNN, in a calculation of the asymmetry. We will see that, although these 
contributions represent almost 30% of the decay width, they affect the asymmetry in a very 

moderate way. 

The work is organized as follows. In Section II we discuss general aspects of the asymme- 
try in the proton emission from the decay of polarized A-hypernuclei. A formal derivation 
of the expressions needed to evaluate the asymmetry parameter is done in Section HI when 
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only lA^ decays are included, resulting in the intrinsic asymmetry a\ , and in Section IV 
when 2A'" decays and FSI are taken into account, resulting in an approximation for the 
observable asymmetry a^. Numerical results are presented and discussed in Section V and 
finally, our conclusions are given in Section VI. 



II. GENERAL CONSIDERATIONS ON THE ASYMMETRY PARAMETER 

Spin-polarization observables for baryon-baryon interactions are important quantities 
which supply additional information to the more usual total cross sections or decay rates 
and thus facilitate the reconstruction of the interaction amplitudes from experimental data. 
For nucleon-nucleon elastic scattering, a complete study of spin-polarization observables 
is given, for instance, in [24]. The formal derivation of the asymmetry parameter for the 
non-mcsonic weak decay of A-hypernuclei is instead provided by [25] . Here we follow a less 
conventional analysis in order to remark some conceptual issues. 

Let us denote with 9 the angle between the momentum of the outgoing proton in the 

— * _ 

Ap — > np weak process and the polarization axis of the hypernucleus. The number of emitted 
protons as a function of 9 can be written as: 

Np{9)^Np^tot{l + Ay{9))/7r, (3) 

where Np^tot is the total number of emitted protons in the decay of the polarized A- 
hypernucleus, while the function Ay{9) introduces an asymmetry in the distribution. By 
construction, it is evident that: 

r d9Np{9)^Np,tot, (4) 

^0 

therefore, 

d9Ay{9)^0. (5) 

Eq. (5) allows one to express Ay{d) as a series of odd powers of cos^. By keeping the first 
term in the series expansion one has: 

Ay{e)^C cos 9. (6) 
This expression is exact for the scattering of two elementary particles, as in the present 

— * 

hadronic description of the Ap — > np weak decay. 
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It is reasonable to write the constant C as the product of the polarization of the hyper- 
nucleus Py times a remaining constant, as follows: 

C = PyAy, (7) 

where the Ay is the hypernuclear asymmetry parameter. Being the A non-mesonic decay in 
a nucleus a complex process, it is evident that also the two-body induced decays Knp — >■ nnp 

— * 

and App — > npp as well as FSI terms contribute to the observable proton number Np{9) of 
Eq. (3). 

If we restrict to the number of protons originated from the elemntary Kp — )■ np process, 
the shell model weak-coupling scheme allows us to make the following replacement: 

PyAy^p^af, (8) 

by introducing the A polarization p\ and the intrinsic A asymmetry a]^ [37]. With these 
definitions, if the weak-coupling limit provides a reliable description of the hypernucleus, 
the intrinsic asymmetry a)^ has the same value for any hypernuclear species. 
We can thus rewrite Eq. (3) as follows: 

^r(^) = Ar;^l(l+PA4^cos^), (9) 

where the index IN refers to the fact that we are considering only the one-nucleon induced 

— * 

decay Ap — > np. Prom this expression, the intrinsic asymmetry is obtained as: 

PA iVpi^(00) + A/"pi^(1800) ■ ^ ' 

Once we consider the two-body induced decay process and FSI as well, the number of 
emitted protons takes the following form: 

Np{e) = Ni^'io) + N^^'+^'^e) . (11) 

As long as Np'^~^^^\9) has a linear dependence on cos 9, it is possible to define an observable 
asymmetry parameter, a]^^'^^'^^^^ , given by a relation which is analogous to the one in 
Eq. (10), which can be compared with the experimental data for the asymmetry a^. 



[37] In the shell model weak-coupling limit it is easy to obtain p\ = — J/ {J + 1) Py for J = Jc — 1/2 and 
PA = Py for J = Jc + 1/2, where J (Jc) is the total angular momentum of the hypernucleus (core 
nucleus). For nuclear matter we have Jc = and then pA = Py. 
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III. FORMAL DERIVATION OF THE INTRINSIC ASYMMETRY 



For computational purposes, we may assume that the hypernucleus is completely polar- 
ized. The intrinsic asymmetry is then given by Eq.(lO) with p\ — Py — 1. 

We now focus on the evaluation of the Np^{9) spectrum. For our practical purpose, we 
can suppose that the polarized A has its spin aligned with the polarization axis (which thus 
coincides with the quantization axis). The evaluation of Np^{6) is rather similar to the 
evaluation of the proton kinetic energy spectrum Np{Tp) described in [20, 22, 23], except for 
two points: i) the angle 9 replaces the proton kinetic energy Tp as variable and ii) we no 
longer sum over the two spin projections of the A, but retain only the up component. This 
second point requires a new evaluation of the spin-summation. To build up an analytical 
expression for Np^{9), let us first express it in terms of the more familiar decay widths as 
follows: 

A^rw=fp(^), (12) 

with Tp{9) = rp(^)/rNM, where Fnm is the total non-mesonic weak decay rate and Fp(^) is 
the proton induced decay rate as a function of 9 [38]. With these definitions, the Np^(9) 
spectrum is normalized per non-mesonic weak decay. 

Before wc give explicit expressions for the ^^-dependent proton spectrum, it is convenient 
to introduce first the weak transition potential: 

yAN^NN^^-^ = J2 aV,^^-^^^(?) , (13) 

T=0,1 

where the isospin dependence is given by 

f 1 for T = 
Or = { (14) 
I Ti • T2 for r = 1 . 

The values and 1 for r refer to the isoscalar and isovector parts of the interactions, 
respectively. The spin and momentum dependence of the weak transition potential is given 
by the function: 

yAN-.NN^^^ = ^G^ml) {Sr{q) <7, q + S'M) ^2 ■ q + PcAq) (15) 

+PL,r(q)(Tl ■q(T2 q + PT,r(g)(0'l X q) ■ ((T2 X q) 

+iSv,ri(l)((^i X <^2) • q} , 



[38] Note that the proton-induced decay rates is obtained as Tp = J^^^^ d9Tp(9) 
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where the index 1 (2) refers to the strong (weak) vertex. The functions Sr{q), S'^{q), Pc,t{q), 
Pl,t{q), Pt,t{q) and Sv,t{q), which include short range correlations, can be adjusted to 
reproduce any weak transition potential. Explicit expressions can be found in [26]. The 
S's {P's) functions are the parity-violating (parity-conserving) contributions of the weak 
transition potential. 

In Fig. 1 we show the Goldstone diagram which has to be evaluated in the calculation of 
Np^{9). The spin summation for this diagram is performed for all particles except for the 




IT>t 

FIG. 1: Direct Goldstone diagram corresponding to the square of the Ap — >■ np decay amplitude. 
A, which is assumed to have spin up. This summation reads: 

sfp\q) = 2 {Sr{q)SAq) + sUq)s'AQ) + PlAq)Pl,Aq) + PcAq)PcAq) (16) 

+2 PtAq)PtAq) + 2 SvA(1)SvAq) 

-2 [Sr{q)Pc,riq) + SAq)PcA<D + S'M)PlA<i) + SA<1)PlA<D 

+2 SvA(i)PTA(i) + 2 SvA(i)PtA(i)\ iA ■ 

It is instructive to note that in the summation over the A spin projection, 

Si:Aiq) + Si:A{q) = 4 {SMSAq) + S'MSAq) + PlA<i)PlAi) + PcAl)Pc,r'{q) 

+2 PTAQ)PT,r'iq) + 2 SvAQ)Sv,r'{q)} , (17) 

the terms between square brackets in Eq. (16) are no longer present. These terms are 
responsible for the asymmetry parameter and arc clearly due to interferences between parity- 
violating and parity-conserving contributions of the weak transition potential in Eq. (15). 

Following [22, 23], we introduce now a partial, isospin-dependent decay width, 
r^^].,(k,kF,9)-> where k is the momentum of the A and kp the Fermi momentum of nu- 
clear matter. This is done for the two isospin channels, r, r' = 0, 1, contributing to the 
spectra. In Fig. 2 we depict the charge-exchange and charge-conserving contributions. The 
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distinction between the two terms is important in the evaluation of A^p^(^) as the kinematics 
of the proton attached to the weak vertex is different from the one outgoing from the strong 
vertex. The partial, isospin-dependent decay widths for the two terms of Fig. 2 are: 




FIG. 2: The two isospin contributions to the direct Goldstone diagram for the Kp — )■ np decay, 
(a) is the charge-exchange contribution, which is not vanishing only for r = r' = 1, while (6) is 
the charge-conserving one. The double arrows represent the A and carry a momentum fc, while 
dashed lines represent the weak transition potential with momentum q. The momentum assigned 
to each intermediate nucleon line is: k — q for the left-most nucleon line {p in (a) or n in (b)), h 
for p~^ and h + q for the right-most nucleon line (n in (a) or p in (b)). 

rg(fc, kp, e) = {Gpmlf^^ j j dqdh S^PHq) 0{qo)0{\k - q\ - k^) (18) 
xe{\h + q\- kF)9{kF - \h\) S{qo - {Er,{h + q) - Er,{h))) 
xS{cos9 - {k — q)z/\k — q\) , 

and 

T%{k,kF,e) = {GFmlfj^^ J J dqdh SSi'^q) 0{qo)e{\k - q\ - kp) (19) 
xe{\h + q\- kF)e{kF - \h\) S{qo - {EN{h + q) - Ej,{h))) 
xS{cos9-{h + q)J\h + q\), 

where the kinematics is explained in Fig. 2. Label (a) refers to the charge-exchange 
contribution (proton attached to the A vertex) and label (b) represents the charge- 
conserving term (proton attached in the strong vertex). In previous equations one has 

qo = — Ej^ik — q) — VV, k^ being the total energy of the A, Ej^ the nucleon total free 
energy and Vn the nucleon binding energy. After performing the isospin summation we 
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obtain: 

r — Ar^"-) \ r^^) \ r^'') r^^) r(^) for\\ 

J- p — 1,1 -r J- 1,1 -r i 0,0 ^0,1 ^1,0) V^'-'J 

where the (fc, kp, ^)-dependence of all functions is omitted for simplicity. Finally, the decay 
rates for a finite hypernucleus are obtained by the local density approximation, i.e., after 
averaging the above partial width over the A momentum distribution in the considered 
hypernucleus, |'0A(fe)P, and over the local Fermi momentum, kpir) — {37r^p(r)/2}-'^/^, p(r) 
being the density profile of the nuclear core. One thus has: 

rp{0) = J dk\Mk)? j dr\Mr)\%{k,kF{r),e) , (21) 

where ■0a(?") is the Fourier transform of ■0A(fc)- The A total energy is given by ko — niA + 
k^/ (2toa) + Va, where Va is a binding potential. 

Finally, by inserting the quantities N^^{e) = Tp{e) = rp(^)/rNM for 0° and 180° in 
Eq. (10) with p\ = 1, the intrinsic asymmetry a]^ is obtained. 



IV. EFFECT OF THE STRONG INTERACTION ON THE ASYMMETRY 

The evaluation of the asymmetry a]^'^'^^'^^^^ , which includes the effects of both 2A'" and 
FSI-induced decay processes, is an involved task and, up to now, analytical expressions were 
given only for the intrinsic asymmetry a]^ , while numerical calculations were performed for 
a]^^^^^ by using the aforementioned hybrid approach incorporating the INC [9, 15]. In this 
section we present for the first time analytical expressions for a]^'^'^^'^^^^ . 

We follow similar steps as in the last section in order to derive Np'^~^^^^{9), which pro- 
vides the total proton spectrum N^^+^^+^^^ (6) = N^^(e) + N^^+^^^{e). This is done by 
introducing the set of Feynman diagrams depicted in Fig. 3 to take care of 2A^ decays and 
FSI effects which result from the action of the nucleon-nucleon strong interaction involving 
the nuclcons produced by the weak decay and nuclcons of the medium. The choice of the 
set of diagrams in Fig. 3 is motivated by previous calculations [20, 22, 23], which show 
that these are the dominant contributions in the evaluation of the nucleon emission spectra. 
Each Feynman diagram is the sum of a number of time-ordering (i.e, Goldstone) diagrams. 
It is in terms of these Goldstone diagrams that one can differentiate among lA^, 2A", pure 
FSI and quantum interference terms (QIT) between lA" or 2A^ and FSI contributions. This 
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FIG. 3: The set of Fcynman diagrams considered in this work for the in-mcdium A self-energy. 
The different time-ordering terms of these diagrams contribute to 2N and FSI-induced decays. 

point is relevant as it shows that, from a quantum-mechanical perspective, each of the above 
processes are included in a unitary description. More details on this point are given in [20]. 

Since in the evaluation with Goldstone diagrams the 2A'" decays are separated contribu- 
tions from FSI-induced decays (which are divided in pure FSI and QIT terms), the 2N+FSI 
proton spectrum reads: 



(22) 



where: 



(23) 
(24) 



Here, f = F/Fnm stands for the decay rate of a particular decay mode normalized per non- 
mesonic weak decay. The functions r„p and Fpp represent the Anp — )> nnp and App — )> npp 
decay processes, respectively, while Fjj represent either pure FSI Goldstone diagrams or 
QIT Goldstone diagrams, accounting for the quantum interference among lA^ or 2N and 
FSI-induced decay processes. The index i in f jj is used to label a particular Goldstone 
diagram obtained from the Feynman diagrams in Fig. 3, while / denotes the final physical 
states of the Goldstone diagram and in the present case can take the values f — np (cut on 
2plh states) and / = npN (cut on 3p2h states), with N = n, p, since we need at least one 
proton in the final state to obtain N^^^'^^'^^^\9). Finally, Nf is the number of protons 
contained in the multinucleon state /. 
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At this point it is necessary to introduce the adopted nucleon-nucleon strong potential: 

y^^(t) = Yl a.^r (t) , (25) 

rjv=0,l 

where t is the momentum carried by the strong interaction, Orj^ is defined in Eq. (14) and 
the spin and momentum dependence of the interaction is given by: 

W = S + ^L,rAt)a, ia,. i (26) 

+VT,r^(t)(<Ti X i) ■ (0-2 X £)} , 

where the functions Vc,tn{^)-i ^l,tn{^) ^^<i '^t,tn{^) adjusted to reproduce any strong 
interaction. 

In the calculation of the diagrams in Fig. 3 the isospin summation is particularly complex 
as one has to differentiate the isospin projection of each particle. We give details on this 
aspect in the present Section and in the Appendix. The main features of the momentum 
dependence of the diagrams were discussed in [20] and references therein. However, the 
important point in the evaluation of the asymmetry is the spin dependence of the diagrams. 

We thus start by considering the spin summation for each Goldstone diagram obtained 
from Fig. 3. This sum is performed for all particles except the A, which again is assumed to 
have spin up with respect to the polarization axis. For the diagrams pp and hh we obtain: 

SZ^Jr'^<l^t) = 2 '^S'^l?) {Vc,r.(t)Vc,.^(t) + V^,.,(t)V^,.^(t) +2VT,r.{t)VT,r'M ' 

(27) 

where 5^*? ^(g) is given in Eq. (16) . The spin summation is more complex for the ph diagram. 
It is convenient to split it in the sum of two terms: 

S'X.r'M^t) ^ Sf^lZr'^rM+Sf^l-^ , (28) 

where 

'55;";;'r(5,t) = 4 {{SMSAq) + PcA<i)PcA<i)) (29) 
+2( Sv,MSvAq) + PT,MPTAq)) yv^^r'^t)} , 
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represents the term which does not contribute to the asymmetry and 

+2( SvAq)PT,Aq) + SvAq)PTAq)) w^^rv (*)> > 

is the term responsible for the asymmetry. In these expressions we have introduced the 
functions: 

y^r^r'. = ^C,r.^Cy^ + ^L,rr,^L,r'^ + 2 Vt, V^, , (31) 

representing the effect of the strong interaction. For simphcity, the t-dcpendence in both the 
W's and V's has been omitted in these expressions. Although Eq. (30) is a more complicated 
expression than the ones in Eqs. (16) and (27), again the asymmetry is originated from the 
interference between parity-violating (5"s) and parity-conserving (P's) terms of the weak 
transition potential. 

The next step is to implement the momentum and isospin summation for each Goldstone 
diagram. In this Section we choose the Goldstone diagram in Fig. 4 as a representative 
example for this evaluation and we leave to the Appendix the remaining contributions. This 
diagram is a particular time-ordering contribution stemming from the hh Feynman diagram 

— * 

in Fig. 3. It contributes to the two-nucleon induced decay mechanisms Anp — > nnp and 

— * 

App — > npp with one and two protons in the final states, respectively. 

Protons in the final state can be in any of the nucleon lines labelled by pi, p2 and p^ in 
Fig. 4. To deal with this matter, it is again convenient to introduce some partial, isospin 
dependent decay widths. For the Anp nnp Goldstone diagram of Fig. 4 we define the 
following rates: 

x^(|fc - g| - kF)e{\h -t\- kp) e{\h -t + q\- kF)e{kF - \h\) 
x9{\h' + t\-kF)e{kF-\h'\) 5{cose-{k-q)J\k-q\) 
xS{qo - {EAh -t + q)- E^(/i) + EN{h' + t) - EAh')) 

S'^^r'-.r^r' Aq,t) 

{EN{h -t)- EN{h) + EN{h' + t)- EN{h')Y 
13 




FIG. 4: Time-ordering term from Fcynman diagram hh in Fig. 3 contributing to 2A'^ decays. The 
momentum carried by the weak transition potential (strong interaction) is q (t). The momentum 
of each particle is k for p\, k — q iov pi, k — t + q iov p2, h + t for P2, h for /12, h' + t for and 
h' for hs- 

^rXlr'. (k, k,,e)^ r;;;f4,,^ (fc, k,, 6) I (33) 

5{cose - (fe - q)J\k - q\) 5 {cos 9 - {h - t + q)z/\h - t + q\) , 



^rXU'. k^,e)^ T'X'.r'. k^, 9) I (34) 

5{cose - (fe - q)z/\k - q\) 5 {cos 9 - {h' + t)z/\h' + t|) , 
each one for the final proton occupying the nucleon line pi, p2 or pa in the diagram, respec- 

— * 

tively. In a similar way, for the reaction Kpp — > npp, we have: 

x^(|fe -q\- kF)e{\h -t\- kp) e{\h -t + q\- kF)e{kF - \h\) 
xe{\h' + t\ - kF)0{kF - \h'\) 

x{5{cose -{k- q)z/\k - q\) + 5{cose - {h' + t)z/\h' + t\))/2 
xS{qo - {EN{h -t + q)- EN{h) + EN{h' + t)- EN{h')) 



{EN{h - t) - EN{h) + EN{h' + 1) - EN{h')Y ' 

rJJ'f.:;:?, (fe, kF. 0) = r'^X'fJ^ i^, kF, e) \ (36) 

5(cos^ - {k - q)J\k - q\) 6 {cos 9 - {h - t + q)z/\h - t + q\) , 
where the sum of the two delta functions in cos^ in Eq.(35) is divided by two in order to 
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retain this multiplicative factor in front of ^pp{0) in Eq. (23). Note that charge conservation 
does not allow particles pi and p2 to be two protons simultaneously. 

^ — * 

The next step is to implement the isospin summation. For the Anp — )■ nnp decay we 
obtain: 

T^hh,pi _ A/rT^hh,pi ,-phh,pi r)-phh,pi\ ('J7\ 
np — 11,11 "•"^ 00,00 11,00 J ' K'^') 

T^hh,p2 p:p^^,P2 I Y^hh,p2 _ 'irvphh,p2 , f-p/i/i,P2 

^ np — oL n il -|- i 00,00 01,11 "I" '-'^ 00,11 

_i_p/i/i,P2 n-phh,p2 n-phh,p2 n-phh,p2 

"T^ 11,00 01,00 11,01 11,01 

r)-phh,p2 I AT^hh,p2 

00,01 "T 01,01 ' 

Y^hh,p3 _ c--p'ift,P3 I -phh.p-^i in-phh.pii ryhhyPi 
^ np — -Jl 11,11 -T- i 00,00 01,11 ^-"^ 00,11 

j_-phh,p3 I nT^hh,p3 _ n-phh,p3 _ n-phh,p3 
"T^ 11,00 "T 01,00 11,01 11,01 

op'i/iiPS /\T^hh,P3 

00,01 01,01 1 

where the (fc, /ci?, ^)-dependence in all the functions are omitted for simplicity. For the 

— * 

Kpp — > npp decay we obtain instead: 

T^hh,pi,p3 _ A(-phh,pi,p3 , -phh,pi,p3 , rfphh,pi,p3\ foQ\ 
PP ~ ^1^11,11 "•"^11,00 "•"^^11,01 ) 1 

T^hh,P2,P3 T^hh,P2,P3 _|_ T^hh,P2,P3 , T^hh,P2,P3 , T^hh,P2,P3 , rfT^hh,P2,P3 

^ PP — ^ 11,11 "T ^ 00,00 "T ^ 00,11 "T ^ 11,00 "T 01,11 

_r)-phh,p2,P3 I r^-phh,P2,P3 , r^-phh,p2,P3 _ /l-nhh,p2,P3 

01,00 "T 11,01 "T ^J- 00,01 01,01 

The last step is to integrate over (fc, kp) in order to implement the local density approx- 
imation as seen in Eq. (21). We have, then: 

r5(^) = v^^^^^ie) + v^^^^-{e) + v^^^^-{e) , (39) 

r^W = rjp^'^i'^'(^^) + rjp^'^^'^^n^)- 

In the Appendix, we show the derivation of some of the other contributions. Once one 
normalizes per non-mesonic weak decay, these expressions are inserted in Eq. (23) and (24) 
to obtain the final result for N^^+^^^{e) [see Eq. (22)]. 

Before presenting the results, we anticipate some elements which emerge from the ob- 
tained analytical expressions and the numerical calculation. First, the ph contribution of 
Fig. 3 turns out to be negligibly small. In addition, we have checked that the behavior of 
jy2N+FSi^Q^ is approximately linear in cos^. We expect this result because for the domi- 
nant hh and pp terms [see Eq. (27)] the spin dependence which generates the asymmetry is 
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given by the same function of Eq. (16) which enters the calculation of the intrinsic 

asymmetry. This allows us to obtain the final expression for the asymmetry as follows: 

Our predictions for Qj^^+2Ar+F5/ -^^ directly compared with the data obtained for the 
observable asymmetry ajy . 



V. NUMERICAL RESULTS 

The weak transition potential y^'^^'^^ of Eq. (13) is described in terms of the usual one- 
meson-exchange (OME), together with the uncorrelated and correlated two-pion-exchange, 
which was shown to have a very important effect on the asymmetry [15]. The OME potential 
is represented by the exchange of tt, rj, K, p, uj and K* mesons within the formulation of [6], 
with values of the coupling constants and cutoff parameters taken from [27] (Nijmegen89) 
and [28] (Nijmegen97f). We present results for both Nijmegen89 and Nijmegen97f weak tran- 
sition potentials for the following reason. The adopted two-pion-exchange potential was in- 
troduced within a chiral unitary approach in [17], together with an important compensatory 
o'-exchange contribution with a a ANou parity-conserving coupling, g%j^^^ — 3.69 G^m^, 
which is the same of the Nijmegen89 potential. At variance, in the Nijmegen97f potential 
one has g%]^^ = 0.17 Gprnl.. One thus expects a difference between the Nijmegen89 and 
Nijmegen97f results for the asymmetry. Although the Nijmegen89-based weak potential was 
the one originally employed in conjuction with the two-pion-exchange mechanism, in a set 
of recent contributions we have used the Nijmegen97f potential and we believe that it is of 
interest to discuss this particular parametrization too. 

For the nucleon-nucleon strong interaction V'^^ of Eq. (25) we have used the Bonn 
potential [29] in the framework of the parametrization of [30] , which contains the exchange 
of n, p, a and u mesons and neglects the rj and 6 mesons. We present results for ]^C, where 
the hyperon is assumed to decay from the lsi/2 orbit of a harmonic oscillator well with 
frequency Huu — 10.8 MeV adjusted to the experimental energy separation between the s 
and p A-levels in [31]. 
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A. Non— mesonic decay rates 

The two-pion-exchange potential is introduced in our microscopic approach for the first 
time here. It is thus important to start our discussion showing the numerical results for 
the non-mesonic weak decay widths. These rates are given in Table I for the two transition 
potentials, Nijmegen89 and Nijmegen97f, without (OME) and with (OME+27r) the two- 
pion-exchangc contribution. Let us start by discussing the independent rates r„, Tp and 
For r„ and r2 all predictions agree with data within error bars; instead, apart from the 
OME result with the Nijmegen97f potential, our predictions overestimate the data for Tp. 
The origin of the agreement for r„ and the disagreement for Fp is not known. However, it is 
the same which leads to a good description of experimental emission spectra involving only 
neutrons and an overestimation of data on spectra involving at least one proton. This is 
proved by the comparison of all the theoretical approaches [19-21] to the single and double- 
coincidence nucleon emission distributions with the corresponding KEK [32] and FINUDA 
[33] data. These nucleon spectra are the real observables in non-mesonic decay, while the 
experimental values of the partial decay rates r„, Fp, etc, are obtained after a deconvolution 
of the FSI effects contained in the measured spectra. The disagreement on the spectra is 
thus the fundamental problem, which also affects the above disagreement on the Fp rate. 

TABLE I: The non-mesonic decay widths predicted for (in units of the free decay rate). The 
most recent data, from KEK-E508 [32] and FINUDA [33, 34], are also given. 



Nijmegen89 Nijmegen97f 



OME OME+27r OME OME+27r KEK-E508 [32] FINUDA [33] FINUDA [34] 



r„ 


0.19 


0.15 


0.16 


0.23 


0.23 ±0.08 


Fp 


0.65 


0.61 


0.47 


0.65 


0.45 ±0.10 


Fi 


0.84 


0.76 


0.63 


0.88 


0.68 ±0.13 


F2 


0.17 


0.26 


0.17 


0.37 


0.27 ±0.13 


Fnm 


1.01 


1.02 


0.80 


1.25 


0.95 ± 0.04 


Fn/Fp 


0.29 


0.25 


0.34 


0.35 


0.51 ±0.13 ±0.05 


F2/FNM 


0.17 


0.26 


0.21 


0.30 


0.29 ±0.13 0.24 ±0.10 0.21 ± 0.07lg;g^ 
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Prom Table I we also see that the effect of the two-pion-exchange potential is different 
when added to the Nijmegen89 and Nijmegen97f OME potentials. This is due to the differ- 
ent values of the (/^Afa; coupling constant previously discussed. While there is a moderate 
reduction of r„ and Tp for Nijmegen89, an increase of these decay rates is observed in the 
case of Nijmegen97f. The behaviour in this later case agrees with what was found in [15, 36] 
with the same weak transition potential [39]. The addition of the two-pion-exchange po- 
tential increases substantially the value of r2 for both the Nijmegen89 and Nijmegen97f, the 
effect being stronger for the later case. There is a certain dispersion among the results ob- 
tained with the different potentials. However, considering the big error bars for data and the 
mentioned discrepancy on the proton emission spectra, we believe that all the four potential 
models of Table I should be also considered in the analysis of the asymmetry parameter. 

B. The asymmetry parameter 

We start by discussing the intrinsic asymmetry a]^ . In Table II we compare our predic- 
tions (first four lines) with the results reported in the literature (last four lines), where the 
updated results of the finite nucleus calculation of [15] have been hsted. 

We obtain a rather sizable asymmetry parameter for the OME Nijmegen89 and Ni- 
jmegen97f models, in agreement with other works, especially with the nuclear matter result 
of [5] . Note that our OME results are more moderate than any of the values found by calcu- 
lations performed in finite nuclei. This is probably due to the more extended Fermi motion 
effects in nuclear matter. The inclusion of the two-pion-exchange mechanism strongly de- 
creases the absolute value of a]^ , especially for the Nijmegen97f model, in agreement with 
what was found by the finite nucleus calculation of Chumillas et a/. [15]. 

Finally, we note that the appreciable difference for the intrinsic asymmetry results eval- 
uated within our two OME -|- 27r models is a consequence of the different values for g'^jvu; 
the Nijmegen89 and Nijmegen97f OME potentials. 

Before moving into the new effects explored in this work, let us comment on the fact 
that our microscopic calculation of the asymmetry parameter takes care of the two isospin 

[39] We note that the results of [15] have been recently revised using more realistic form factors [36] . Although 
the numerical values have slightly changed, the qualitative aspects of the two-pion-exchange mechanism 
remain the same. 
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TABLE II: Theoretical determinations of the intrinsic asymmetry parameter for ]^ C. The calcu- 
lations reported from the literature were performed within shell model approaches, except for the 
nuclear matter result of Dubach et al. The resuls of the model of Chumillas et al. [15], correspond 
to the ones updated in [36] . 



Model 




OME (Nijmegen89) 


-0.39 


OME (Nijmegen89) + 27r 


-0.23 


OME (Nijmegen97f) 


-0.35 


OME (Nijmegen97f ) + 27r 


-0.071 


Ref. and Model 


Dubach et al. [5], OME (NM) 


-0.44 


Parreno and Ramos [7], OME 


-0.55 to -0.73 


Barbero et al. [8], OME 


-0.53 


Chumillas et al. [15, 36], OME 


-0.48 


Chumillas et al. [15, 36], OMEH 


-27r -0.0062 



channels depicted in Fig. 2. In fact, these contributions are automatically encoded within 
the antisymmetric character of the final two-nucleon wave-function used in finite-nucleus 
calculations of the weak decay. However, the diagrammatic approach employed here is useful 
in the sense that it allows one to keep track of the importance of the different contributions 
to the asymmetry parameter. It turns out that diagram (a) is the dominant contribution. 
This can be explained as follows. Let us denote with p„ (p^) the momentum carried by the 
proton in diagram (a) {{b)). Prom the kinematics of diagrams (a) and (6) we have: 

Pa^k-q^-q, (41) 
Pb = h + q^q, 

where, for the purpose of this explanation, it is a good approximation to assume that k and 
h are much smaller than q. The different sign but similar magnitude of and means 
that a negative asymmetry from the charge-exchange diagram (a) is reduced in magnitude 
by a positive asymmetry from diagram (6). The competition between the diagrams (a) and 
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TABLE III: Effect of the nucleon-nucleon strong interaction on the asymmetry parameter for j{C. 
In addition to the contribution from the Ap — )■ np weak decay leading to the intrinsic asymmetry, 
we consider the action of 2N and FSI-induced decays. The most recent data are also shown. 







Nijmef 


?en89 


Nijmegen97f 


Eth (MeV) 


Asymmetry 


OME 


OME+27r 


UME 


(JME+ZTT 





"a 


—0.386 


—0.225 


-0.352 


-0.071 





1N+2N 

"a 

IN+FSI 

"a 

1N+2N+FSI 

"a 


-0.366 
-0.184 
—0.234 


-0.212 
-0.009 
—0.071 


-0.318 
—0.043 
-0.132 


-0.063 
0.082 
0.032 


30 


1N+2N 

(Z A 

^N+FSI 

"a 

^1N+2N+FSI 

"a 


-0.355 
-0.149 
-0.196 


-0.197 
-0.003 
-0.056 


-0.307 
-0.034 
-0.115 


-0.060 
0.096 
0.037 


50 


1N+2N 

IN+FSI 
"A 

„1N+2N+FSI 
"A 


-0.319 
-0.123 
-0.156 


-0.149 
0.019 

-0.018 


-0.255 
0.014 

-0.058 


-0.049 
0.112 

0.069 


KEK-E508[35] 








-0.16±0.28lg:oo 



(6) produces a reduction in the absolute value of the asymmetry. This type of analysis will 
be particularly useful for the 2N and FSI effects discussed below. 

In Table III we present our predictions for the asymmetry when 2N and FSI-induced 
decays are considered together with IN decays. Since any experiment is affected by a kinetic 
energy threshold for proton detection, Eth, results are also given for different values of Eth- 

The values obtained for the asymmetry depend on two important effects. The first one 
is the dynamics of the weak transition. In particular, one can consider or not the two- 
pion-exchange potential. This has been analyzed in detail in [9, 15] , and our results confirm 
those findings. Moreover, the asymmetry depends on what we call "kinematic effect" . The 
introduction of the different 2A'" and FSI contributions enlarges the available phase, leading 
to a particular kinematics for each contribution; the weight imposed by the nucleon-nucleon 
strong interaction on the different kinematics (and also the restrictions due to Eth) niodifies 
the relation between Np{0^) and A'p(180°). The microscopic model is particularly suitable 
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for the study of this kinematic effect. Note that the division between the dynamic and 
the kinematic effects is possible because the spin summation representing the interference 
between parity-violating and parity-conserving terms of the transition potential has the 
same expression, given by Eq. (16), for the intrinsic asymmetry and for the dominant 2N 
and FSI contributions to the observable asymmetry. 

It is instructive to recall that the value of the asymmetry parameter is a consequence of 
a delicate balance between parity-conserving and parity-violating amplitudes governed by 
the dynamics of the weak decay mechanism, and also depends on the phase space allowed 
for the emitted nucleons which might be enhanced or decreased in some places by strong 
interaction effects or kinematical cuts. Any new contribution to the decay process will intro- 
duce changes in the number of protons emitted parallel, iVp(0°), and antiparallel, iVp(180''), 
to the polarization axis, therefore affecting the value of the asymmetry which is determined 
by the difference Np{Q^) - A/p(180°) measured relative to the sum A^p(0°) + A^p(180°), as seen 
in Eq. (40). It is therefore illustrative to represent the function Np{9) as a function of cos^, 
as seen in Fig. 5 for the OME Nijmegen89 model, including the IN induced decays (dot- 
ted line), adding the 2A'"-induced modes (dashed line), adding only FSI effects (dash-dotted 
line), and finally incorporating all the contributions together (solid line). Similar plots are 
obtained for the other three potential models employed in this work. 

It is clear that the 2A'" mechanism enhances the number of emitted protons at all emitted 
angles but having a slight preference for directions opposite to the polarization axis, hence 
the size of the slope of the dashed line in Fig. 5 is a little bit larger than that of the dotted 
line. This would increase the magnitude of the asymmetry but the larger number of protons 
gives finally rise to a slight decrease, as seen in Table III. Diagrammatically speaking, 
we know that, for lA^ decays, the asymmetry receives the main (negative) contribution 
when the proton is attached to the A vertex (diagram (a) of Fig. 2). The other (positive) 
contribution, with a neutron outgoing from the A vertex (diagram (6) of Fig. 2), tends to 
reduce the absolute value of the asymmetry. In the case of 2A" decay diagrams, one has 
npp and nnp final states, where the proton(s) can be located at the A vertex or in any of 
the two remaining positions. It is the increased number of positions for the final proton(s) 
that produces a further reduction (although small) of the asymmetry parameter. We note 
that the small effect ol2N decays on the asymmetry parameter corroborates the assumption 
done in [9] . 
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0.4 




FIG. 5: N^'^ (dotted line), N^^+^^ (dashed line), N^^+^'^^ (dash-dotted line) and total 
j\j-iN+2N+FSi (^gQ[[f^ line), as functions of cos0, in the the case of the OME Nijmegen89 weak 
transition potential. 

As far as FSI effects are concerned, we observe in Fig. 5 that they remove antiparallel 
protons and, on the other hand, more strength is added at parallel kinematics. Since the 
total number of protons is almost unchanged, this reduction of slope observed for the dot- 
dashed line also translates in a subtantial decrease in the magnitude of the asymmetry. In 
order to analyze further the different FSI contributions, it is convenient to write: 

N^^'{9) = N^P'\e) + N^P^'^ie) , (42) 

where each of the two terms on the rhs receive contributions from each of the diagrams in 
Fig. 3, by cutting on 2plh or 3p2h states, respectively. By construction, N'^^^^{9) originates 
from a QIT between lA^ and FSI-induced decays, while Np''^^{9) may come either from a 
27V-FSI QIT term or from a pure FSI-induced decay. The microscopic model allows us to 
inspect the behavior of each term. We find that the term Np^'^^{6) is positive-definite and 
has a similar behavior to the one already discussed for Np^{9), i.e. slightly more protons 
are emitted antiparallel to the polarization axis. On the contrary, N^^^^iO) turns out to 
be negative while its kinematic behavior is very similar to the lA^-induced decays, which 
produce a large negative asymmetry parameter. Therefore, the effect of the negative NpP^^{9) 
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contributions goes in the direction of inverting this behavior, giving rise to a subtantial 
decrease in the size of the asymmetry or even reverting its sign, as in the case of the 
Nijmegen97f+27r model. 

We now pay attention to the behavior of the asymmetries of Table 111 with the enegy 
cut Eth- We observe that the size of the asymmetry a]^'^'^'^ decreases slightly for increas- 
ing Eth- This reduction can be explained from a microscopic point of view by inspecting 
the momentum distribution predicted by our approach for the three particles, pi (particle 
outgoing from the A vertex), p2 and ps, stemming from 2A'" decays, shown in Fig. 12 of [21]. 
Particles are named in that figure with the same notation as in Fig. 4 of the present work. 
The distributions for particles pi and ^3 are very similar to each other and are peaked at a 
lower momentum than the distribution for p2. Due to isospin reasons, the main (negative) 
contribution to the asymmetry is obtained when a proton is located in pi, while protons 
in P2 and/or p^ reduce the magnitude of the asymmetry. The effect of Eth is to reduce 
the importance of the particle pi with respect to p2- This explains the reduction in the 
magnitude of a]^~^'^^ for increasing Eth- Note, however, that the decrease is much stronger 
in the case of the a]^~^^^^ asymmetry, and this is also the behavior of the complete cal- 
culation, a]^'^^'^'^^^^ . In order to understand this behavior, we recall that the energy cut 
removes nucleons, making Np{0^) + A^p(180°) smaller and, consequently, the magnitude of 
the asymmetry larger. But the final consequence of the Efh cut on the asymmetry will be 
determined by whether the increase in size due to the removal of nucleons is counterbal- 
anced by the changes in the slope A/p(0°) — A^p(180°). In Fig. 6 we show the effect of this 
cut for N^^+^^{9) (dashed lines) and N^^+'^^+^^\9) (solid lines) as functions of cos ^, for 
the OME Nijmegen89 model. We clearly see a reduction in the number of protons as well 
as a decrease in the slope with increasing Eth ■ These two effects modify the asymmetry 
parameter in opposite ways and the results of Table III show that, within our models, the 
reduction of the asymmetry due to the decrease in the slope dominates over the increase 
associated to the removal of particles. The reduction in the slope is much more pronounced 
for the FSI contributions. The final result is that we observe a substantial reduction in the 
magnitude of the asymmetry a]^'^'^^'^^^^ with an increasing energy cut. 

This behavior contrasts with the INC results of [9, 15], where, for increasing Eth, the 
size of the asymmetry a]^~^^^^ increases and tends to the intrinsic value. We note that 
the INC model for FSI originates from a semi-classical description which has an intuitive 
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IN 

1N+2N 




FIG. 6: The functions N^^'^'^^ (dashed hnes) and jy^^+'^^+P'Sl (gohd hnes) for different energy 
cuts, Eth = 0,30 and 50 MeV, in the the case of the OME Nijmegen89 weak transition potential. 
The dotted line corresponds to N^'^ at Efh = MeV. 

interpretation. Nucleons are tracked in their way out of the nucleus as classical particles. 
Sometimes a nucleon leaves the nucleus without any interaction with the medium, in other 
cases it scatters one or more times with the other bound nucleons. Therefore, a nucleon 
emerging from an elementary non-mesonic decay can change momentum, direction and 
charge, other nucleons can be emitted as well, etc. Clearly, the random character of these 
FSI processes is responsible for the strong reduction of a]^~^^^^ by about a factor two or 
more with respect to the intrinsic asymmetry a]^ [9, 15]. The introduction of an energy cut 
Effi affects mainly those nucleons which have suffered scattering processes. For increasing 
Eth, the nucleons coming from elementary decays (not affected by FSI) become dominant 
and the asymmetry tends to the intrinsic vahic. This is also reflected in Fig. 2 of [9] by the 
tendency of Np^~^^^\9) 1^^^ to move towards Np'^ {9) as the energy cut is increased, both 
functions becoming very similar (in size and slope) at around Eth — 50 MeV. 

The situation for our microscopic approach is different as it is based on quantum mechan- 
ics, where QIT play an important role. It was shown in [20] that the 2plh and 3p2h terms 
of the proton kinetic energy spectra N^^^Tp) = N^P^^iJp) + N^P'^''{Tp) have a different 
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behavior from each other. While NpP'^^{Tp) gives a positive distribution, has its maximum 
for Tp — and decreases for increasing Tp, the QIT N'^P^^{Tp) is a negative bell-shaped 
distribution with the minimum at Tp = 80 MeV. Thus, a non-vanishing energy cut E^h ap- 
preciably reduces N^^"^^ while leaving A^^^^'* almost unchanged, and this later contribution 
is the one producing a significant decrease in the slope of N^^^'^^^^^\9). Therefore, as Et^ 
is increased, the magnitude of the asymmetry parameter decreases. 

We end our discussion by comparing our results with experiment. The asymmetry data 
reported in Table III was obtained by KEK-E508 for a kinetic energy threshold Eth of 
about 30 MeV. Despite the noticeable differences among the whole set of predictions, they 
arc all compatible with experiment due to the large error bar of data. By considering our 
results for Eth = 30 MeV and the central value of the experimental data, the best agreement 
is obtained for both Nijmegen89 and Nijmegen97f models when the two-pion-exchange 
potential is not included. Certainly, our calculation shows that the effect of the two-pion- 
exchange is very important in asymmetry calculations, but to establish definite conclusions 
on the effect of this potential more detailed studies are required. In general, the addition 
of any new contribution to the weak transition potential has to be done consistently with 
the rest of the potential itself (which might require some readjustment to reproduce the 
observables) and with the approach adopted in the calculation. However, in order to obtain 
fruitful information from these studies new and more precise data is needed to constrain the 
unknown parameters of the weak decay models. 



VI. REMARKS AND CONCLUSIONS 



We have discussed a microscopic diagrammatic formalism to evaluate the asymmetry in 
the distribution of protons emitted in the non-mesonic decay of polarized hypernuclei. The 
calculation is performed in nuclear matter and then extended to finite hypernuclei (a^C) by 
means of the local density approximation. Our approach takes into account both the 2A'" 
decay mechanism and the nucleon FSI in a unified many-body scheme. The effect of the 
2A'" decays on the asymmetry parameter is evaluated here for the first time. The present 
work is also the first one to implement the FSI on the asymmetry parameter by means of a 
quantum-mechanical microscopic approach. In addition to the usual OME weak transition 
potentials, which we take from the Nijmegen89 and Nijmegen97f parametrizations, we have 
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also considered the effect of the two-pion-exchange potential introduced in [17]. We give 
results for both the intrinsic asymmetry parameter, a]^ , and for the asymmetry parameter 
modified by the 2A^-induced mechanisms and FSI effects, a]^'^^^'^^^^ , which is the one that 
can be compared to the observed asymmetry a^. 

While the effect of 2N is predicted to be rather limited, the nucleon FSI turned out to 
be very important: they reduce the magnitude of the asymmetry parameter, making all the 
weak transition potential models adopted in this work capable of describing consistently 
the experimental data for and for the non-mesonic weak decay rates. In particular, the 
large error bars in the observable asymmetry do not allow us to determine which of the two 
potential models, OME or OME+27r, provides the best description of the experiments. 

To the best of our knowledge, the only former work which evaluated the intrinsic asym- 
metry in nuclear matter is due to Dubach et al. [5], where an approximate scheme rather 
different from ours (neglecting 2A'" decays and nucleon FSI) was employed. The action of 
FSI was considered within a semi-classical description in [9, 15], by means of an INC model. 
In a former calculation [20] , it was shown that the INC model and the present microscopic 
approach provide similar results for the nucleon emission spectra in the non-mesonic weak 
decay of unpolarized A hypernuclei. Our results for a]^ and a]^^^'^^ with a vanishing pro- 
ton kinetic energy cut, Efh — 0, fairly agree with each other too. However, the situation 
changes for non-vanishing values of Efh- For increasing Ef^, the negative asymmetry of the 
INC model increases in magnitude, while a decrease is observed in the microscopic model. 
One should note that the two schemes represent rather different approaches to the prob- 
lem of dealing with nuclear correlations after the weak decay takes place. The microscopic 
model of the present work provides a reliable method that can be improved systematically. 
It however ignores multinucleon processes that are accounted for, semiclassically and via 
multi-step processes, in the INC model. To determine which is the most realistic approach, 
an accurate experimental determination of the asymmetry parameter, possibly exploring its 
E'tft-dependence, would certainly be welcome. 

One should always keep in mind that the main motivation in the study of the non-mesonic 
weak decay of hypernuclei is to extract information on strangeness-changing baryon-baryon 
interactions. The understanding of the F^/Fp ratio and the asymmetry parameter suggests 
that a fairly reasonable knowledge of non-mesonic decay has been achieved. However, we 
have obtained agreement with all the experimental data employing different parametriza- 
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tions of the weak transition potential. Due to the lack of precise data for the asymmetry 
parameter, we find that the role of the two-pion-exchange mechanism, which was essential 
to reproduce this observable in some models [15], can not even be firmly established here. 
In any case, what is certain is the agreement with the set of data can only be achieved after 
a proper development of approaches that take care of nucleon FSI. Due to the special nature 
of the in-medium non-mesonic weak decay, these are complex models, but they are required 
to establish a link between theory and experiment. 

Finally, we recall that there still remains an important disagreement between theory an 
experiment for the hypernuclear non-mesonic weak decay: theoretical evaluations of nucleon 
emission spectra involving protons strongly overestimate the experimental distributions. 
This discrepancy may not be isolated but hidden behind the errors bars in the data for the 
decay rates. An additional aspect that has not yet been studied but which could lead to a 
non-neghgible contribution to the nucleon spectra is the inclusion of the A(1232)-resonance 
in our many-body Feynman diagram scheme. We intend to study this problem in the future. 
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Appendix 

Here we present the explicit expressions needed in the evaluation of N^^ {0) starting from 
the Feynman diagrams pp and ph in Fig. 3. We omit the derivation of the expressions for 
Np^^{9) obtained from the same diagrams, as they can be obtained from the N^^ [9) ones 
after some simple changes: the spin-isospin structures are the same, as well as the general 
expressions, except for some step functions and energy denominators. 

We begin with the contribution of the diagram pp of Fig. 3. First, we introduce the 
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partial, isospin-dependent decay widths for Rnp — > nnp: 

rf;,- ^ (g^_L_^ j,„ j,„ I „n jih'OM (43) 

x^(|fc -q\- kF)e{kF -\h + q\) e{\h -t + q\- kF)e{kF - \h\) 
xe{\h' + t| - Uf) e{kF - \h'\) S{cose - (fc - q)J\k - q\) 
x5{qo - {EN{h -t + q)- EnUi) + Ejv(/i' + t) - EN{h')) 

{EM{h -t + q)- E^ih + q) + E^ih' + t)- E^ih'))^ ' 

rf;a.v(fe' kF, e) = r™ .,^(fe, kF, e)\ (44) 

5{cose -{k- q)z/\k - q\) 5(cos^ - {h - t + q)z/\h - t + q\) , 

rf;.f4.v(fe' kF, 0) = r™ ...(fe, kF, e)\ (45) 

5(cose - (fc - g) J|A; - q\) S{cose - {h' + t)J\h' + t|) , 
where pi, p2 and pa indicate the position of the final proton. In a similar way, for the reaction 

— * 

Kpp — > npp we have: 

x^(|fe -q\- kF)e{kF - \h + q\)e{\h-t + q\- kF)e{kF - \h\) 
x9{\h' + t\-kF)e{kF-\h'\) 

x(5(cose- {k - q)z/\k - q\) + 5(cos^ - {h + q - t)J\h + q- t\))/2 
xd{qo - {EN{h -t + q)- E^ih) + E^ih' + t) - E^ih')) 

{E^ih + q-t)- Ej^ih + q) + E^ih' + t)- E^ih'))^ ' 

rf;{rrv kF,9) = rf;,p/,^^ {k, kF, e) \ (47) 

5{cose -{h + q- t)z/\h + q-t\)^ 5 {cos 9 - {h' + t)J\h' + t\), 

rf;^f^/.?. {k, kF,e) = rf;,p/,^^ {k, kF, e) \ (48) 

5(cos^ - (fe - q)z/\k - q\) 5(cose- {h' + t)J\h' + t|) . 
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The next step is to implement the isospin-summation to obtain: 

np — ^y'- 11,11 ^ 00,00 ^ 11,01 ) 1 V^^) 
^ np — JJ- 11,11 -r i 00,00 "T ^ 11,00 "T 00,11 

_OpPP,P2 I nyPP,P2 _ r)'pPP:P2 _ r)j^PP:P2 , AT^PP,P2 

11,01 't' 01,11 00,01 01,00 "T 01,01 ' 

^ np — 11 11 -t- i 00,00 "T ^ 11,00 "T 00,11 

_9PPP,P3 _ f;pPP,P3 _ OpPP,P3 I OpPP,P3 _ y|pPP,P3 

11,01 01,11 00,01 "T 01,00 01,01 ' 

where the (fc, kp, ^)-dependence of all these functions has been omitted for simplicity. In a 

— * 

similar way, for the Kpp — > npp reaction we have: 

rPP,Pi,P2 ^ i6r™'^% (50) 

PPP,P1,P3 _ A(-pPP,Pl,P-A I pPP,Pl,P3 _ 9pPP,Pl,P3\ 
PP ~ ^l-*- 11,11 "I" -"^ 11,00 ^-"^ 11,01 )l 

ppp,P2,P3 _ pPP)P2,P3 I pPP,P2,P3 I pPP,P2,P3 , 9pPP,P2,P3 _ 9pPP,P2,P3 
PP ~ ■'-11,11 00,00 T" 00,11 11,01 ^■'-01,11 

_1_9PPP,P2,P3 9pPP,P2,P3 /1pPP,P2,P3 

"T^^ 00,01 01,00 01,01 

The final point is to employ Eq. (21) to implement the local density approximation. We 
have, then: 

r^^(^) = ' (^) + (^) + ^np"' {0) , (51) 

rpp(e) = vp^^^^'P^e) + r^^'^'i'P3(^) + rpp'P^'P3(^) . 

Finally, the pp contribution to Np^{9) is obtained by Eq. (23). 

We then consider the 2A'" decay contribution from the ph diagram of Fig. 3. We follow 
the same steps of the former contributions. We start by introducing the partial, isospin- 
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dependent decay widths for Anp — > nnp: 

x^(|fc -q\- kF)e{\h -t\- kp) e{\h -t + q\- kF)0{kF - \h\) 

xe{kF -\h + q\) e{\h' + t\- kp) e{kF - \h'\) S{cose - (fc - q)J\k - q\) 

xS{qo - {EN{h -t + q)- EN{h) + Enih' + t) - EN{h')) 

^ ^lT';Tr,T'M^'^) 

EN{h - t) - EN{h) + EN{h' + t) - EN{h') 

1 

X- 



EN{h + q-t)- EN{h + q) + E^ih' + t)- EN{h') 

rS;f4,,^(fe, kF. 0) = KX\^,jk, kF, e)\ (53) 

5{cose - (fe - q)J\k - q\) 5 {cos 9 - {h - t + q)J\h - t + q\) , 

r5f4rv(fe' kF, 0) = rS;f4,,^(fe, kF. 9)\ (54) 

5{cose - (fe - q)J\k - q\) 5 {cos 9 - {h' + t)J\h' + t\) , 
where pi, p2 and pa indicate the position of the emitted proton. In a similar way, for the 

— * 

reaction App — > npp we have: 

x^(|fc - q| - - t| - kF) e{\h -t + q\- kF)e{kF - \h\) 

xe{\h' + t\-kF)e{kF-\h'\) 

x{5{cose - (fe - q)z/\k - q\) + 5{cose - {h' + t)J\h' + t\))/2 
xd{qo - {EN{h -t + q)- Ej^{h) + EN{h' + t)- EN{h')) 

X " 



{EN{h -t)- EN{h) + EN{h' + t)- EN{h'))^ ' 



rSf.^^?. (fc, kF. 0) = KX^r'. (k, kF, 9) I (56) 

5{cose- {k - q)z/\k - q\) 5 {cos 9 - {h - t + q)z/\h - t + q\) , 
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The next step is to implement the isospin summation to obtain: 



np 




T^ph,p2 

np 


— "J^ 11,11 

"r"-'- 01,11 


pp/l,P3 

np 


_ qpP^, P3 
— '^^ 11,11 

01,11 ■ 



"I" ^ 11,00 "I" 11,01 J > ) 
T" ^ 00,00 ' 00,11 "T ^ 11,00 

_ C)j^pll,P2 _ 9pp/t,P2 _ r)pp/l,p2 _|_ /|ppA,P2 

01,00 11,01 ^-"^ 00,01 "T 01,01 ' 

I pp/iiPa I cpp/iiPa I pp/i,P3 
"T ^ 00,00 "T 00,11 "T ^ 11,00 

I OpP'liP3 _ 9pP/l,P3 _ r)pP'j,P3 _ y|pP/l,P3 

"T 01,00 11,01 00,01 01,01 ' 

where the (fc, kp, ^)-dependence of all these functions has been omitted for simplicity. In a 

— * 

similar way, for the Kpp — > npp reaction we have: 

ppfe,pl,P3 _ /|/'_pp/l,Pl,P3 I pP'j,Pl,P3 _ 9pp/l,Pl,P3\ /roN 
PP ~ ^ 11,11 "T ^ 11,00 11,01 ) ' V^°/' 

ppfe,P2,P3 _ pP'*:P2,P3 I pP^,P2,P3 I pP'l,P2,P3 , pP'*,P2,P3 _ 9pPft,P2,P3 
PP ~ 11>11 00,00 "T -"^ 00,11 "T 11,00 01,11 

_9pp/l,P2,P3 I nyph,P2,p-i I 9pP'l,P2,P3 _ /lpP^,P2,P3 

01,00 11,01 -r 00,01 01,01 ' 

One thus has to perform the local density approximation, through Eq. (21), to obtain: 

viio) = r^'^H^) + + rt^'^'id) , (59) 

and finally the ph contribution to N!^^(9) is obtained by Eq. (23). 
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